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algorithm are compared with the CRB and the comparison leads to validate this estimator. 



Abstract 

This paper presents a least mean square (LMS) algorithm for the joint estimation of acoustic and mean flow velocities from 
laser doppler velocimetry (LDV) measurements. The usual algorithms used for measuring with LDV purely acoustic velocity or 
mean flow velocity may not be used when the acoustic field is disturbed by a mean flow component. The LMS-based algorithm 

CL, 

allows accurate estimations of both acoustic and mean flow velocities. The Cramer-Rao bound (CRB) of the associated problem 
is determined. The variance of the estimators of both acoustic and mean flow velocities is also given. Simulation results of this 
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Laser Doppler Velocimetry for Joint Measurements 
of Acoustic and Mean Flow Velocities : LMS -based 
Algorithm and CRB Calculation 

I. Introduction 

Laser Doppler Velocimeter (LDV) is an optical technique allowing direct measurement of local and instantaneous fluid 
velocity. This method is nonintrusive and is based on optical interferometry for estimating the velocity of scatterers suspended 
in a fluid by means of the frequency analysis of the light scattered by the seeding particles [1]. 

For fluid mechanics measurements, the particle velocity can be considered as constant during the transit time of the seeding 
particle through the measurement volume (defined by the interferometry fringes volume) and the frequency of the LDV signal 
is constant during this period [16]. Typical order of magnitude of mean flow velocities are from a few meters per second up 
to higher than the acoustic celerity (supersonic flow). The data processing consists then to estimate the power spectral density 
(PSD) of the velocity signal, from Poisson-based randomly distributed samples. PSD may be estimated by interpolating the 
randomly distributed samples, by resampling the interpolating signal and by compensating the effect of interpolation in the 
Fourier domain [5]- [20]. The autocorrelation function (ACF) may also be reconstructed from the randomly distributed samples 
and the Fourier transform of the estimated ACF gives an estimation of the PSD [11]. Lastly, Kalman filtering may be used for 
estimating the PSD [2]. 

For sine acoustic excitation, the particle velocity is no longer constant and the LDV signal is frequency modulated [21]- 
[10]. To estimate the particle velocity from these signals, specific signal processing techniques are used as spectral analysis [6], 
[22], [25], photon correlation [18] or frequency demodulation associated to post-processing methods [8], [23], [24]. Typical 
order of magnitude of mean flow velocities are from a few micrometers per second up to 100 millimeters per second, for 
frequencies in [10 - 4000] Hz. 

On one hand, for most acoustic measurements, the particle velocity can be considered as the sum of an AC-component due 
to acoustic excitation and a weak DC-contribution due to flow. When the particle oscillates in the measurement volume during 
further acoustic periods, the effect of the flow can be reduced and usual post-processing methods may be used [24]- [9]. 

On the other hand, the DC-flow component prevents in many cases the use of the post-processing methods given by [8], 
[23], [24], because the signal time length is less or largely less than one acoustic period. The aim of this paper is to estimate 
both the dc (flow) and ac (acoustic) components from such LDV signals. 

Lazreq and Ville [13] measured the acoustic velocity in presence of mean flow by means of a probe consisting in a hot wire 
and a microphone. Their results showed a good agreement between the theory and the experiment but this probe cannot be 
considered as nonintrusive. LDV has also been used by adapting the slotting technique to estimate the acoustic particle velocity 
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in a turbulent flow [15] with a 2D-LDV velocimeter. The acoustic impedance was estimated by means of a LDV probe and 
with a microphone probe and the different results were compared. Finally, Boucheron et al [4] has developed a new method 
of signal processing called 'perio-correlation' in order to estimate sine acoustic velocity in strong mean flow by LDV. 

In this work, the sine acoustic excitation is supposed to be perfectly known and a frequency demodulation technique [8] 
is performed to estimate the particle velocity from the LDV signal. In this paper, we propose a new method to estimate 
jointly the acoustic particle velocity (amplitude and phase) and the mean flow velocity from the velocity signal. This method 
is based on the least mean square (LMS) algorithm. The mean flow velocity, the amplitude and phase of acoustic particle 
velocity are estimated for each seeding particle crossing the measurement volume. Furthermore, the Cramer-Rao bound (CRB) 
of the associated problem is calculated. The CRB gives the lowest variance of any unbiased estimator and consequently yields 
theoretically the minimum uncertainties linked to the velocity estimations (acoustic and mean flow velocities). Lastly, simulated 
data are processed, in order to validate the LMS-based algorithm and to compare the variance of the results with the Cramer-Rao 
bound. 

Section II deals with the LDV principles including the velocity signal modeling and the associated signal processing for 
acoustic applications. In section III, the data processing based on the least mean square algorithm is explained and the Cramer- 
Rao bound of both the mean flow and acoustic velocities are determined. Finally, the results of the Monte Carlo simulation 
are shown and compared to the Cramer-Rao bounds in section IV, for acoustic frequencies in [125 — 4000] Hz, for acoustic 
velocities in [0.05 — 50] mm.s -1 and for mean flow velocities in [0.05 — 5000] mm.s -1 . 

II. Fundamentals of Laser Doppler Velocimetry 

In this section, we consider time-varying signals such that t G [t q — T q /2, t q + T q /2], t q being the central time of the signal, 
T q being a time of flight, and q being associated to a given seeding particle. 

A. Laser Doppler Velocimetry Principle 

In the differential mode, two coherent laser beams are crossed and focused to generate an ellipsoidal probe volume, in which 
the electromagnetic interferences lead to apparent dark and bright fringes [1]. 

The velocity v q (t) of the seeding particle denoted q is related to the scattered optical field due to the Doppler effect. The 
light intensity scattered by the particle crossing the probe volume is modulated in amplitude and frequency. The frequency of 
modulation F q (t) is called Doppler frequency and is given by 

^) = ^ = ^ S in(0/2), (1) 

where v q (t) is the velocity of the particle along the x-axis, i the fringe-spacing expressed as a function of the angle between 
the incoming laser beams and their optical wavelength Al (Fig. Q]). 

The diffused light is collected by a receiving optics and is converted into an electrical signal by a photomultiplier (PM). 
This signal can then be modelled as [8] 



s q (t) = A q (t)(M + cos <£,(*)), 



(2) 
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Fig. 1. Optical setup of LDV system. When the particle q crosses the measurement volume, the light is scattered in all directions and the burst signal s q (t) 
is collected by the photo detector. Data processing of s q (t) allows then to estimate the mean flow and particle acoustic velocity. 

where M takes into account the positive sign of Cramer-Rao Bound (CRB) the light intensity. In (gj, the amplitude modulation 
linked to the normally distributed light intensity across the beam section is written as 

A q (t) = K 9 e-V ,d *W (3) 

where K q is related to the laser beam, the PM sensitivity, the electronic amplification, the observation direction and the 
scattering efficiency of tracer q. Furthermore, (3 is related to the probe geometry and d q (t) is the projection of the time-varying 
particle displacement along the a;-axis in the probe volume. Similarly, the phase modulation in © is described by 

9 ( t ) = 2^^+^ > (4) 
i 

where <fio is the initial phase due to optical setup. Furthermore, we denote x q (t) the signal such that 

x q (t) =s q {t)+w(t), (5) 

where w(t) is the additive noise [24]. 

In order to avoid any ambiguity on the sign of the velocity, a Bragg cell tuned to frequency Fb = 40 MHz is used to shift 
the frequency of one of the lasers. The signal s q (t) is consequently written as 

s q (t) = A q (t)(M + coa{2nF B t + 2nd q {t)/i + cp )). (6) 

The offset component M is then canceled by an high-pass filtering and the signal s q (t) is down shifted to zero thanks to a 
quadrature demodulation (QD) technique [14]. The actual signal, called burst signal, can finally be written as 



s q (t) = A q (t) cos(27rd q (t)/i + <j> ). (7) 

B. Doppler Signal Modeling in Acoustics 

Considering only pure sine acoustic waves and supposing that the mean flow velocity is constant inside the probe volume, 
the projection along the a>axis of the velocity of a particle q subjected jointly to the sine acoustic wave and the mean flow 
field can be expressed as 

Vq(t) = V c . q + V ac COs(2lTF ac t + OC ), (8) 



where v CtQ is the mean flow velocity of particle q, V ac and <fi ac are the amplitude and phase of the acoustic particle velocity and 
F ac is the known frequency of the pure sine acoustic excitation. The amplitude modulation of the burst signal (O associated 
to the particle q may be written as 

A q (t) = K q exp[P{v c . q (t - t q ) + — ^- sin(27rF„ c i + <M)] 2 - (9) 

Similarly, the phase modulation (01 of the burst signal associated to the particle q is 

<f> q (t) = —v c , q {t - t q ) + sm{2irF ac t + OC ). (10) 

i 2nr ac 

We note that the flow velocity v c , q can change from a particle q to another while the acoustic parameters v ac and <^> ac 
are independent of q. Thus, when the acoustic wave is disturbed by a mean flow, assuming that the particles q cross the 
measurement volume at different random central times t q , without time overlapping between bursts q and q + 1, the Doppler 
signal can be written as 

8(t)=^2a g (t)=A D (t)ooa[4> D (t)], (11) 
where the amplitude and phase respectively express as 



A D {t) = < 



A q (t), te[t q -T q /2,t q + T q /2] 
0, otherwise, 



(12) 



and 



4> q {t), te[t q ~T q /2,t q + T q /2] 
){t) = \ (13) 

0, otherwise. 



Furthermore, the time of flight of the tracer q is defined as [7] 

V2D X 

T q = * *, (14) 

Vc.q 

where D x is the length of the probe volume in the x-axis and the associated number of acoustic periods is 

V2D X 

AW = 'Fac- (15) 

Vc,q 

As expected, the fastest the particle crosses the probe volume, the lowest the time of flight and the number of acoustic periods. 
An example of a typical Doppler signal is shown on Fig. 12a), where the different particle times of flight are associated to 
different mean flow velocities. 

C. Doppler signal processing 

The aim of the signal processing developed after sampling the Doppler signal is to estimate jointly and burst-by-burst the 
acoustic particle velocity (amplitude V ac and phase <f> ac ) and the mean flow velocity v c i9 . This procedure is usually split into 
two stages. After a detection procedure [7], a frequency demodulation of the Doppler signal s(t) is performed by using a 
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Fig. 2. (a) Example of a Doppler signal, (b) Associated velocity signal. Burst 1 is associated with a low mean flow velocity corresponding to N per acoustic 
periods largely higher than 1. Burst 2 is associated with a high mean flow velocity corresponding to Nper acoustic period largely lower than 1. Burst 3 is 
associated to a mean flow velocity corresponding to N per < 1 acoustic period. 

time-frequency transform to estimate the instantaneous frequency F q (t), or equivalently ([TJ the velocity signal v q (t), burst by 
burst [23]. Note that the detector selects only bursts corresponding to one tracer in the measurement volume. Secondly, the data 
processing of the estimated velocity signal v q (t) allows to obtain both components of the acoustic and mean flow velocities 
for each burst. This first stage is described in this subsection and the second one (data LMS-based processing) is explained in 
the section HVl 

According to ([T), the velocity signal associated to the particle q expresses as 



v q (t)=iF q (t), 



(16) 



where i is the fringe spacing. Thus, the problem consists to estimate the mean value v c<q , the amplitude V ac and the phase <f> ac 
of the estimated velocity signal associated to each burst q from the actual noisy burst signal x q (t). Fig. [2] shows an example 
of an noiseless simulated Doppler signal (a) and the associated velocity signal for three non-overlapping bursts (b). 

In section [TTTJ the Cramer Rao bound of the problem is calculated. Then, a method based on a least mean square algorithm 
is presented in section iPVl and is applied to simulated velocity signals v q {t) in section |V1 



III. CRB CALCULATION 

We recall that the Cramer-Rao bound (CRB) gives the lowest bound of the variance an unbiased estimator may reach (if 
it exists) [12]. As explained in [12], the CRB alerts us to the physical impossibility of finding an unbiased estimator whose 



variance is less than the bound. In the case of single tone signals, CRB were calculated by Rife and Boorstyn [17] in 1974. 
The CRB of LDV signals were also studied in the case of fluid mechanics [3]- [19]. In the case of sine acoustic excitation, 
the CRB of LDV signal were also studied by Le Duff [14]. 

We focus here on the problem of calculating the Cramer-Rao bound (CRB) of the following problem. The velocity data are 
assumed to be such that 

u[n] = v[n; 9] + w[n], (17) 
for n € [no, rii], where w[n] is the WGN, w[n) ~ A/"(0,ct 2 ), the data being modeled according to 



v[n; 9} = v c + V ac cos(27r f ac n + cj> ac ), 



(18) 



where f ac = F ac /F s , F s being the sampling frequency, v c = v c _ q , and where the unknown parameters are gathered in 



= \v B V r , 



c v ac V^acl 



(19) 



We furthermore suppose that f ac ^ and f ac ^ ~. 

A. Cramer-Rao bound (CRB) for one burst 

The CRB is given by the inverse of the Fisher information matrix 3(9), CRB(#) = 3(9) , where the Fisher information 
matrix is given by [12] 



3(9) 



kl 



dv[n; 9) dv[n; 9} 



cr 2 ^ d9 k d9i 

n—no 



for k. I € [1, 3], for 6 — [v c V ac 0ac] T - The derivatives in (|20] >, according to (fT8l >, lead to 



3(9) 



V 



cos(/3)sin(7Af) v ae sin( / 9)sin(7jV) 

sin(7) sin(7) 

cos( / 3)sin(7jy) jv , cos(2j3)sin(27Af) v ac sin(2 / 3)sin(27Af) 

Sin(7) 2 + 2Sin(27) 2Sin(27) 

v ac sin(/3)sin(7jv) v ae sin(2^)sin(27jv) Nv£ r _ v a 2 c cos(2 ) 3)sin(27Af) 

sin(7) 2sin(27) 2 2sin(27) 



where AT = m — no + 1, and where 



7 = Kfac 



(3 = 2irf ac n + irf ac (N - 1) + (j> a 



(20) 



(21) 



(22) 
(23) 



We define the linear signal-to-noise ratio (SNR) of the velocity signal as 

V 2 

SNR=-t, 



(24) 
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and we then have upon inversion 



N 2 - 



sin(2 7 /y) \ 



Vac V Sm(27) / 

4 SNR 

2 2 I sin(27) y I sin(7) / 1 sin(27) \ sin(7) 



var( Vc ) > CRB(« C ) - -± -{ — , (25) 

N3 N ( sin(2 7 AT) \ N ( sin( 7 jv) \ , sin(2 7 AQ ( sin( 7 jv) 



iV 2 - Ar cos(2/3 )5i!M _ 2sin 2 (/3) sirurjv) 
V 2 v A ; sin(2 7 ) KH > \ sin( 7 ) 



2 



var(F ac ) > CRB(K C ) = - 2 - 2 * f (26) 

A N s N I sin(2 7 jy) \ N I sin( 7 jV) \ srn(2 7 AQ / sin( 7 /Y) 



2 2 1 sin(2 7 ) / I sin( 7 ) J sin(2 7 ) I sin( 7 ) 

/Y 2 + jVcos(2/3)^« _ 2cos 2 (/3) f siru^iv) 
I w ' sin(2 7 ) \ sin 7 ) 
> CRB((f> ac ) = ttt^ ; ; — r ri • ( 27 > 



2 SNR v 

N 3 
2 



n ( sin(2 7 Af) A N f sin( 7 jy) A sin(2 7 AQ ( sin( 7 AQ 
2 I sin(2 7 ) J I sin( 7 ) y sin(2 7 ) I sin( 7 ) 



5. Cramer-Rao Bound (CRB) for Nb bursts 

We now assume that the algorithm developed in MI- Al is used for estimating the unknown parameters 9 = [v c V ac 4>ac] T , in 
the case of Nb bursts. The main difference between this problem and the one developed above is that the index no is not anymore 
a constant, but might be modeled as a discrete random variable, uniformly distributed in [0, N ac ], where N ac = nint F s /F ac , 
nint[] being the nearest integer. As a consequence, the discrete random variable (3 given in ( |23T l, which appears in ( |26ll27l ) 
is uniformly distributed in [tt(N — l)/ ac + <Pac, 7r(/Y — l)/ ac + 4>ac + 27r]. Averaging the terms linked to j3 in ( |26] i and ( |27l ) 
consequently leads to 

< cos(2/3) >=< sin(2/3) >= 0, (28) 

and 

< cos 2 (/3) >=< sin 2 (/3) >= i. (29) 

This finally yields 

2 

sin( 7 /Y) 



yy2 _ 

V 2 \ sin( 7 ) 
var(K c ) > CRB(K C ) = 7-^- - 2 - (30) 

z jya jy / sin(2 7 /v) \ ^ / sin( 7 AQ \ sin(2 7 jV) / sin( 7 JV) 
sin(2 7 ) J \ sin( 7 ) y 1 sin(2 7 ) ^ sin( 7 ) 



and 2 

sin( 7 AQ 



1 V Sln (7) 

2 SNR 



var(0 ac ) > CRB(0 ac ) = ; ^ ^ ^ -j. (31) 



n± _ n( sin(2 7 jv) V _ sin( T jy) V sin(2 7 iv) / sin( 7 Af) \ 
2 2 I sin(2 7 ) J I sin( 7 ) y sin(2 7 ) I sin( 7 ) y 

In the following, we use the expressions (l25l l for w c and ( f30b for 14 c for studying the CRB of the problem. We recall that 
N depends on v c (l36l l. As a consequence, the CRB of v c d25l > and the CRB of V ac <f30b both depend on v c and 14 c , while the 
CRB of <fi ac (l3TT > is independent of V ac . 



C. Asymptotic behavior of Cramer- Rao Bound (CRB) 

In Appendix ITU we give the expressions of the asymptotic CRB of 9, for both cases 2^N <C 1 (N per <C 1/(2tt)) and 
2 1 N > 1 (N per > 1/(2tt)). 

In the asymptotic case 2^N <C 1, we prove (1631164b that the relative variance of v c and V ac are 

var(t; e ) CRB(v c ) 1 45 1 v 3 c V* c 

vl ~ vl SNR^ 2 7/2^5 Fs FL > ^ > 

var(Kc) CRB(Kc) _ 1 45 1 v 5 c 

Both CRBs of v c and V ac are proportional to v\V2 c and inversely proportional to F^ c . Consequently, doubling the mean 
flow velocity yields an 15 dB increase of the variance of both v c and V ac . Similarly, doubling the amplitude of the acoustic 
particle velocity V ac leads to a 6 dB increase of the variance of both v c and V ac . Lastly, doubling the frequency of the pure 
sine acoustic wave leads to a 12 dB decrease of the variance of both v c and V ac . We also note that doubling the length of the 
probe volume D x yields a 15 dB decrease of the variance of both v c and V ac - 

In the asymptotic case 2^N 3> 1, we prove that (l72H73b 

^ (34) 

and 

var(Kc) > 77^ 7=r— v c Vl c . (35) 

V acj _ ^ Dx p s ac 

Thanks to the exact d25l [30l [3Tb and asymptotic ( I32II35I ) expressions of the CRB, the minimum uncertainties linked to the 
velocity estimations (acoustic and mean flow velocities) are completely known. In section (lIVI l. the LMS-based algorithm is 
introduced. It is then applied in section ([V) to simulated data in order to be compared with the CRB. 

IV. Least mean square algorithm 

From a practical point of view, the actual velocity signal is uniformly sampled. Consequently, the number of samples N q 
associated to the particle q is derived from ( fl4l i. as 

V2D X F S 

= X ^ 5 (36) 

V c ,q 

and the associated number of acoustic periods ( fTBI ) is now defined as 

V2D X 

N pe r = ~ ~F ac . (37) 

V c ,q 



The sine-wave fit is then solved by minimizing the cost function V{6), 

n— no 

with respect to the unknown parameters 6 < [T9K where u[n] and u[n;0] are respectively given by < fT7] > and ( fTSI t, and where 
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N = nx — no + 1. In the Appendix 1, the equations d50ll52l ) respectively give the expression of v c , a ac = V r ac cos(0 ac ) and 



b ac = V ac sin((/> ac ) as a function of u and f ac . Once a ac and b ac are estimated, the unknown acoustical parameters of 9 express 
as 



■ , b 2 

ac 1 ac' 



atan-^ 



(39) 



V. Numerical results and discussion 

In this section, we compare the CRB with the LMS-based algorithm developed in|lV] According to the values of the acoustic 
and mean flow velocities to be analyzed, the following values for F ac and V ac are chosen : 



F ac G [125 250 500 1000 2000 4000] Hz, 
V ac S [0.05 1.58 50] mmr 1 . 



(40) 
(41) 



The phase <p ac is supposed to be equal to ir/4, and we use an adimensional parameter a v for the value of v c , such that 



= — G [0.01 0.05 0.1 0.5 1]. 



(42) 



For each numerical simulation, the sampling frequency is F s = 350 kHz, the probe volume length along the x— axis is 
D x = 0.1 mm and 10000 bursts are analyzed. The simulator is performed by Matlab. 

Fig. [3] to Fig. [5] show typical results of the relative variances \ai(v c )/v 2 (a) and var(V ac )/V r a 2 c (b), for the different values 
of F ac with comparison to the theoretical CRB of v c ( l25l ) and of V ac i26i . Each figure is related to a given value of F ac . 
Furthermore, for each value of F ac , three sets of signals are analyzed, each set corresponding to one of the bursts of Fig. [2] 
respectively N per ^> 1 (Burst 1), N per > 1 (Burst 2) and N per < 1 (Burst 3), where N per is the number of acoustic periods, 

N per = Nfac- 
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Fig. 3. Comparison of the relative variances of v c (a) and Vac (b) estimated by a LMS algorithm (continuous) with the theoretical CRB (dashed), for 
F ac = 125 Hz. Bursts (1 - 3) refer to Fig.f2] (Burst 1) : V ac = 1.58 mm.s -1 , a v = 0.1, v c = 15.8 rnm.s" 1 . (Burst 2) : V ac = 50 mm.s _1 , a v = 0.1, 
v c = 500 mnu" . (Burst 3) : Vac = 50 mm.s — , a v = 1, v a = 50 mm.s -1 . 



First of all, the LMS-based estimator is near the theoretical CRB, so that we can maintain that this estimator is efficient. 
Moreover, the relative variance of v c is weaker than the one of V ac (except for a v — 1). Indeed, for values of v c and V ac such 
that a v = Vac/vc < l/y/2 as CRB(w c ) ~ CRB(V r ac )/2, we have 
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Fig. 4. Comparison of the relative variances of v c (a) and V ac (b) estimated by a LMS algorithm (continuous) with the theoretical CRB (dashed), for 



F ac = 500 Hz. Bursts (1 - 3) refer to Fig. |2] (Burst 1) : V ac = 50 mm.s" 
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Fig. 5. Comparison of the relative variances of v c (a) and V ac (b) estimated by a LMS algorithm (continuous) with the theoretical CRB (dashed), for 
Fac = 4000 Hz. Bursts (1 - 3) refer to Fig. [2] (Burst 1) : V ac = 50 mm.s" 1 , a v = 0.1, v c = 500 mm.s -1 . (Burst 2) : Vac = 50 mm.s" 1 , a v = 0.01, 
v c = 5000 mm.s -1 . (Burst 3) : Vac = 50 mm.s -1 , a = 0.05, v c = 1000 mm.s -1 . 



CRB(iQ CRB(Kc) 



ac 



Moreover, the values of the relative variances of v c and V ac drastically depend on the values of v c , V ac and F ac as shown 
by (f25) and (|3D). 

- For N per > 1 (Burst 1), the relative variances of v c and V ac are respectively in [—70, —42] dB and [—52, —30] dB. 
Such estimations may consequently be considered as very accurate. 

- For Np er "C 1 (Burst 2), the relative variances of v c and V ac are respectively in [—27, 1] dB and [10, 20] dB. The estimation 
of v c is accurate enough for low values of v% V^ c /F^ c (l32l . while the estimation of V ac is clearly unacceptable, whatever 
the parameters v c , V ac and F ac d33l l. 

- For AT per < 1 (Burst 3), the relative variances of v c and V ac are respectively in [—70, —30] dB and [—30, —23] dB. Such 
estimations may also be considered as very accurate. 

Furthermore, for velocity signals with time length largely lower than one acoustic period, we can use the asymptotic case 
expression of CRB of v c and V ac . Giving a maximum value of relative error, respectively E Va for v c and E Vac for V ac , we 
consequently have 

CRB(v c ) 1 45 1 vlV?; c ^ 

■ < E Vc (44) 

J ' x 1 - s ± ac 
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and 

CRB(Kc) = 1 45 1_± <F (A r, 

VI SNR 7r 4 29/ 2 D X F S F± c ~ Vac ' 1 ; 

As a consequence, for a given set of setup known parameters D x , F s and F ac , we may give the maximum values v% V Q 2 C 
and have to reach for yielding an error less than respectively E Vc and E Vac . 

Lastly, from the expressions of the CRB of v c d25l l and V ac (l30l l, we can calculate the number of acoustic periods the time 
length of the velocity signals may have, for leading to an error less than a given value E. Tables HlHIIl give a summary of such 
results. Each table corresponds to a given burst of Fig. [2] 

TABLE I 

Number of acoustic periods N per {for» c and for V ac ) leading to an error less than E(%) for F ac = 125 Hz and V ac = 50 mm.s -1 . 



SNR (dB) 


10 


20 


30 


E (%) 


0.1 


1 


10 


0.1 


1 


10 


0.1 


1 


10 


Nper (V c ) 


> 0.8 


> 0.3 


> 0.2 


> 0.45 


> 0.2 


> 0.09 


> 0.3 


> 0.12 


> 0.05 


N per (Vac) 


> 10 


> 0.75 


> 0.25 


> 5 


> 0.5 


> 0.2 


> 0.8 


> 0.25 


> 0.1 



TABLE II 

Number of acoustic periods N per (for-u c and for V ac ) leading to an error less than E(%) for F ac = 500 Hz and V ac = 1-58 mm.s -1 . 



SNR (dB) 


10 


20 


30 


E (%) 


0.1 


1 


10 


0.1 


1 


10 


0.1 


1 


10 


Nper («c) 


> 0.4 


> 0.16 


> 0.06 


> 0.25 


> 0.1 


> 0.04 


> 0.18 


> 0.06 


> 0.02 


Nper (V ac ) 


> 10 


> 6 


> 0.4 


> 10 


> 0.75 


> 0.25 


> 5 


> 0.9 


> 0.15 



TABLE III 

Number of acoustic periods N per (foru c and for V ac ) leading to an errorless than E(%) for F ac = 4000 Hz and V ac = 50 mm.s -1 . 



SNR (dB) 


10 


20 


30 


E (%) 


0.1 


1 


10 


0.1 


1 


10 


0.1 


1 


10 


Nper (V c ) 


> 0.6 


> 0.25 


> 0.1 


> 0.4 


> 0.15 


> 0.05 


> 0.25 


> 0.1 


> 0.03 


Nper (V ac ) 


> 20 


> 5 


> 0.6 


> 10 


> 4.5 


> 0.04 


> 10 


> 0.7 


> 0.25 



For example, table|T]may be read as follows. To obtain a relative error for v c less than 0.1 % for SNR= 10 dB, the minimum 
number of acoustic period for the velocity signal is 0.8. In the same way, to obtain a relative error for V ac less than 1 % for 
SNR= 20 dB, the minimum number of acoustic period for the velocity signal is 0.5. Tables ITTHTin give the minimum number 
of acoustic periods the velocity signal should have for obtaining relative errors less than 0.1 %, 1 % and 10 %, for SNR equals 
to 10 dB, 20 dB and 30 dB for 500 Hz and 4000 Hz respectively. 

As expected, the mean flow velocity v c is estimated with a great accuracy from a very low number of acoustic period. For 
example, to obtain a relative error of 1 % for v c , the number of acoustic period is always less than 0.3 whatever the SNR, F ac 
and V ac - On contrary, the results for the estimation of the acoustic velocity are much more contrasted. For a SNR of 30 dB, 
the estimation of V ac associated with a relative error less than 1 % is possible for a number of acoustic period N per > 0.9. 
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But, when the SNR is less than 30 dB, the number of acoustic periods associated with a relative error less than 1 % may be 
largely bigger than 1. 

The tables also show the influence of the acoustic frequency on the estimation of the particle acoustic velocity. The higher 
the frequency, the higher the number of acoustic periods for an accurate estimation of V ac . For a relative error equals to 10 
%, the number of acoustic periods are the same whatever the frequency. But for a relative error equals to 1 % or 0.1 % the 
estimation of the particle acoustic velocity is easier for a low frequency. On contrary, the influence of the frequency on the 
estimation of the mean flow velocity is the opposite. The higher the frequency, the lower the number of the acoustic periods 
for an accurate estimation. 

With regard to the results of these study, a three-steps new approach can be proposed to improve the estimation of the 
acoustic particle velocity in presence of mean flow. The first step consists in the estimation of the mean flow velocity for each 
burst with the least mean square (LMS) algorithm. Then, the estimation of the mean flow velocity may be subtracted from 
the velocity signal. Finally, a "rotating machinery" technique associated with a synchronous detection allows to estimate the 
acoustic particle velocity with a great accuracy [9]. 

VI. Conclusion 

A new method for estimating jointly the acoustic particle and the mean flow velocities from a LDV signal is presented. It 
is based on the least mean square (LMS) algorithm and it performs well in the estimation of the velocities. The performance 
of the method has been investigated by means of numerical tests and the results of the simulation have been compared to the 
Cramer-Rao bounds of the associated problem. It is shown that the LMS-based estimator is near the theoretical CRB, so that 
the estimator is efficient. 
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Appendix I 
Derivation of LMS problem 



Inserting ( fl8l l into d38l l leads to 



1 711 ( \ 2 
V{6) = — f u n - (w c + a QC cos(27r/ ac n) + 6 ac sin(27r/ QC n) J , 



(46) 



where 



a ac = V ac cos((f> ac ), 
b ac = V ac sin(</> oc ). 



(47) 



Solving the following linear problem 



9V(8) _ n 
dV ac u ' 



allows to write analytically the unknown parameters. In the following, we note 



(48) 



n-ni n=n\ ^ / n=n\ 

D = — ^ cos 2 (27r/ ac n) ^ sin 2 (27r/ ac n) - — ( ^ cos(27r/ ac n)sin(27r/ ac n) 



n— no 
n-ni 



n— no 
, n— ni 



—3 ^2 sin 2 (27r/ ac ?i)[ ^2 cos(2vr/ QC n) ) --3 ^ cos 2 (27r/ ac n) I ^ sin(27r/ ac n) 



n—rtQ ^ n— no n— no 

n-ni n-ni n=m 



+ jV3 cos ( 27r /ac^)sin(27r/ ac n) ^ cos(27r/ oc n) ^ sin(27r/ ac n). 



(49) 



The mean flow velocity i? may then be written as 



_^ / n—rii n-ni n— ri\ n=m \ n-ni 

3n( ^2 sin(27r/ ac n) ^ cos(27r/ ac n)sin(27r/ ac n) - ^ cos(27r/ ac n) ^ sin 2 (27r/ ac n) j ^ u„cos(27r/, 



iV 3 L> 



+ lv^n( 51 cos ( 27r /ac") 51 cos ( 27r /a C n-)sin(27r/ ac n) - ^ sin(27r/ QC n) ^ cos 2 (27r/ oc n) ) ^ u n sin(27r/ a 



n— no 
n— ni 



iV 3 D 

1 

7V 3 L> 



n— no 
n— ni 



n— no 
n— ni 



n— no 
n— ni 



n— no 
n— ni 



n— no 
n— ni 



n— no 
n— ni 



^ cos 2 (27r/ ac n) ^ sin 2 (27r/ QC n) - ( ^ cos(2tt f ac n)sm(2ir f ac n) \ j ^ 



(50) 



Similarly, the acoustic parameters express as 



7V3£)( X! ^ 2 {2irf ac n) f sin(27r/ ac n) J w n cos(27r/ ac n) 

^ n=no n— no ' ' n— no 

, n— ni n— ni n— ni 



n— ni 



2 v n= 



n— ni 



3— I ^ cos(27r/ ac n) ^2 sin(27r/ QC n) - ^ cos(27r/ QC n)sin(27r/ ac n) ) ^ u„sin(27r/ Q 



iV 3 L> 
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^ y n—ri\ n—ni n-ni v n— ni 

3n( 51 sin(27r/ oc n) cos(27r/ ac n)sin(27r/ ac n) - cos(27r/ ac n)sin 2 (27r/ QC n) J u„, (51) 



N 3 D 



n—riQ 



and 



&ac = ^J^f X! cos2 ( 27r /ac"-)f cos(27r/ ac n) J J ^ u n sm(2irf ac n) 

^ n— no 7i— no n— no 

^ y n— ni n— ni n— ni v n— ni 

+ JPD \ ^ cos ( 27r / Qc71 ) ^2 sin ( 2lr /aci) - X! cos ( 27r /ac"-)sin(27r/ QC n) j w„cos(27r/ ac n) 

* n— no n— no n— no n— no 

y n— ni n— ni n— m v n-ni 

+ ( cos(27r/ ac n) ^ cos(27r/ ac n)sin(27r/ ac n) — ^ sin(27r/ ac ?'i)cos 2 (2^/ ac n) ^ u„. (52) 

^ n— no n— no n—no n—no 

Appendix II 
Asymtotic CRB 

In this Appendix, we write the CRB of v c (|25T >. V ac (f30b and ac (|3~TT > respectively in both asymptotic cases 

2 1 N < 1, (53) 

and 

2 7 A > 1, (54) 

where 7 and N = N q are given by (l22l and (f36b . Using (f36l > and (l22l . we note that d53l and (l54i i are respectively equivalent 
to 

2V2irD x F ac < u c , (55) 

and 

2V2nD x F ac > u c . (56) 

Firstly, we suppose that 27A <C 1 and that 7 <C 1 which means that the actual velocity signal corresponds to largely less 
than one acoustic period. The Taylor expansion at the 7th order of the sine functions in d25l ), ( f3Qb and ( f3Tb respectively yields 

45 1 

var^)^-^--, (57) 
var(K c )>^^, (58) 

var(0 oc ) > ^ c ^k c W- (59) 
Using ( 136b and d22l . we note that (157l|59t may respectively be written as 

45 1 v 5 

M^^^pr^, (60) 



2 , ti— r, 
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Writing (|24) into ( I60II62I I leads to 



45 1 v 5 

vaKK^-^ — (61) 
45 1 v 5 

-(^)>- 2 ^^pr^%- (62) 



1 45 1 ullV? 
var(<; c ) > — __^£, (63) 



SNR 7T 4 2 7 /2 




p4 ' 

ac 


1 45 


1 


v 5 V 2 

u c ¥ ac 


" SNR tt 4 2 9 / 2 


•■D*F e 


F A 

ac 


1 45 

> 


1 





SNR tt 4 2 9 / 2 D\F e Ft 



(64) 
(65) 



where SNR is the linear signal-to-noise ratio. 

Secondly, we now suppose that 27^ ^> 1 which means that the actual velocity signal corresponds to largely great than one 
acoustic period. The asymptotic CRB is then such that 



var( t ; c ) > — , (66) 



N 

and 



2cr 2 

var(F ac ) > — (67) 



2 



2(7 

var(0 ac ) > —. (68) 



Using (l36*T l. we note that d66*1l6*8l may respectively be written as 



D x F e 
and 



var(t;c) - uwj:- m 

var(Kc) > ^rv c (70) 



var(<M > (71) 

L 'i r e v ac 



Lastly, inserting d24l i into ( 1691711 ) finally leads to 



var(w c ) > — ^J- v c V 2 , (72) 

V c> ~ SNR 2 3 / 2 D x F e QC ' 

var(V ac ) > — — ^ v c V 2 (73) 

V ac) _ ^jy^ ac 



16 



and 

var(0 ac ) > ^^JTVc (74) 
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